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Abstract: The dynamics of finite temperature U(N) gauge theories on S 3 can be described, 
at weak coupling, by an effective unitary matrix model. Here we present an exact solution to 
these models, for any value of N, in terms of a sum over representations. Taking the large 
N limit of this solution provides a new perspective on the deconfinement transition which is 
supposed to be dual to the Hawking-Page transition. The large N phase transition manifests 
itself here in a manner similar to the Douglas-Kazakov phase transition in 2d Yang-Mills 
theory. We carry out a complete analysis of the saddle representation in the simplest case 
involving only the order parameter Till. We find that the saddle points corresponding to 
thermal AdS, the small black hole and the large black hole can all be described in terms of 
free fermions. They all admit a simple phase space description a la the BPS geometries of 
Lin, Lunin and Maldacena. 
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1. Introduction 

Though we know of many instances where a gauge theory is a holographic description of a 
gravitational theory, we are yet to understand the precise way in which a local diffeomorphism 
invariant theory in one higher dimension is encoded in the dynamics of the gauge theory. In 
a sense, the redundancy of diffeomorphisms has been largely eliminated in the gauge theory 
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description. But this has come at the cost of losing information about the locality of the bulk 
description. Is there a natural way in the gauge theory to restore the redundancies which 
characterise the geometrical description of the bulk? 

A partial hint comes from the beautiful work of Lin, Lunin and Maldacena [1] who showed 
that the geometry of a class of half-BPS solutions of the bulk theory is completely fixed by 
specifying a single function of two of the bulk coordinates. This function, which takes values 
either zero or one in the entire two dimensional plane, was identified with the phase space 
distribution of free fermions describing the half BPS dynamics in the gauge theory. In other 
words, the configuration space of the bulk, with its redundancies, was identified with the 
phase space of the boundary degrees of freedom. In fact, the quantisation of the space of 
BPS configurations on the gravity side agrees with those of the free fermions [2] [3]. Notice 
that the fermionic phase space description is also a redundant one since it is the shape of the 
perimeter of the "filled fermi" sea that completely determines everything. The phase space 
picture therefore appears to be a step in the right direction. 

However, the half-BPS case seems to be very special and the picture of free fermions is 
not likely to be generally applicable. It is therefore a bit of a surprise that, in this paper, 
we find a similar free fermion phase space description in the non-supersymmetric context of 
finite temperature AdS/CFT. It is very well known [4] that the thermal partition function of 
the gauge theory exhibits a behaviour which is qualitatively similar to the Hawking-Page [5] 
phase diagram on the gravity side. In particular, we find a free fermionic description in the 
weakly coupled gauge theory, for each of the saddle points that correspond to Thermal AdS, 
the (unstable) "small" AdS Schwarschild black hole as well as the "big" AdS black hole. In 
each case there is a simple region in phase space which is the filled fermi sea. 

Our starting point is the effective unitary matrix model that describes the holonomies of 
the Polyakov loop at weak gauge coupling [6] [8]. As was argued in [8], in the free U(N) Yang- 
Mills theory on S 3 at finite temperature, all modes are massive and can be exactly integrated 
out, except for the zero mode of Aq. The dynamics of this interacting mode is naturally 
expressed in terms of a unitary matrix model for the holonomy U, along the thermal S 1 . At 
weak coupling, we can continue to integrate out all the other modes and end up with (a more 
complicated) effective matrix model for U. These matrix models have been well analysed, 
in the large N limit, in terms of the collective field a (6) which is the eigenvalue density of 
U. They have been shown to exhibit a phase structure which describes the deconfinement 
transition and is qualitatively very similar to that of the Hawking-Page description of AdS 
gravity at finite temperature [8] [9] [10]. 

In this paper, we present an exact solution to the partition function of these matrix 
models, which is valid for any finite N. The answer is in terms of characters of the conjugacy 
classes of the symmetric group with a sum over different representations and classes. While 
explicit, the expressions are, in general, quite complicated. At large N we expect the answer 
to show the non-analytic behaviour, as one varies the temperature, which is characteristic 
of a phase transition. This is seen in our expressions from the fact that at large N, there 
is a dominant saddle point in the sum over representations. The nature of this saddle point 
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exhibits non-analytic jumps as one varies the temperature. This is similar to how the large 
N phase transition of Douglas-Kazakov [12], in 2d Yang-Mills theory, manifests itself. 

The quantitative method of analysis, as in 2d Yang-Mills, introduces a density u(h) for 
the Young tableaux that label the representations. This essentially measures the number of 
boxes in the rows of the tableaux. One can write an effective action for u(h) in the large 
N limit and study its saddle points. We do this analysis very explicitly for the simplest 
and physically important case 1 where one has only terms involving TrU and TrC/L The 
actual saddle point equations are close to that of models studied by Kazakov, Staudacher and 
Wynter [13], though those cases did not exhibit a phase transition. In our case, one finds, not 
surprisingly, exactly the phase diagram obtained by the usual eigenvalue density analysis. 

However, what is of interest in the present analysis, is the nature of the saddle point 
representations uo(h), in both the low and high temperature phases. It turns out that they 
bear a simple relation to the saddle point eigenvalue densities ctq{9). Essentially the two 
turn out to be functional inverses of each other. The best way, in fact, to state the relation 
between the two is to view (h, 9) as coordinates on a two dimensional phase space and define 
an appropriate region R with constant fermion density p(h, 9) = ^ in its interior and zero 
outside. It then turns out that 



In other words, the region R is determined by the shape h (9) which is obtained from inverting 
the equation 9 = iru$(h). 

It turns out that the region R corresponding to thermal AdS is given by the unit disk in 
phase space. This is indeed what one also obtains in the LLM picture [1] for the global AdS 
spacetime. The regions corresponding to the small black hole and and the big black hole are 
more complicated kidney-shaped geometries as shown in Figs. 7-9. The shape of these regions 
is not modified in functional form when one includes perturbative corrections in terms of an 
effective action involving only the relevant winding number one modes. Thus at least in the 
weak coupling expansion this geometry of the phase space distribution is robust and therefore 
can be expected to capture some essential features of the corresponding bulk geometries. It 
would be very interesting to learn what these features might be. In particular, it is natural 
to ask whether there is a direct translation into a supergravity solution like in the LLM case. 
Gaining an understanding of these points might help us learn why the matrix models capture 
the dynamics of the gravity phase transition so well. 

The plan of the paper is as follows. In the next section (Sec. 2) we review the unitary 
matrix models that describe the finite temperature dynamics at weak coupling. We also 
recapitulate the results that follow from a large N analysis in terms of the eigenvalue density 
a (9) and the correspondence with the phase diagram on the gravity side. In Sec. 3 we write 

x The (a, b) model studied in [10] falls, for instance, in this class. 




(1.1) 
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down the exact finite N solution to the models at zero coupling and also show how the method 
of solution generalises to the weakly coupled case. In Sec. 4 we analyse the large N limit of 
the exact solution for models involving only TrU and TrU^ . We do this in terms of the Young 
Tableaux density u(h) and find a phase transition as expected. We find the expressions for 
the saddle points uo(h) that dominate at both low and high temperature and compute their 
free energies to find agreement with the results of the eigenvalue density analysis. In Sec. 5.1, 
we show how the results of Sec.4 imply a relation of uq(K) with the saddle point eigenvalue 
densities (Jq(9). We show how this relation can be simply understood in terms of a free 
fermionic phase space picture in Sec. 5. 2. In Sec. 6 we close with various comments on the 
possible implications of these results which need to be fleshed out in the future. Appendices 
contain details of some of the calculations as well as some generalisations. 

2. The Finite Temperature Partition Function for Gauge Theories on S 3 
2.1 The Effective Action for the Holonomy 

In the AdS/CFT correspondence the four dimensional gauge theory lives on the boundary 
S 3 (together with the R direction for time) of the AdS§ spacetime. In the finite temperature 
version, the (Euclidean) gauge theory now has a thermal S 1 instead of R. Studying the 
dynamics of the thermal gauge theory on S 3 x S 1 offers some important simplifications. In 
the free gauge theory (defined as the A = 9ym^ ~^ limit), most of the modes are massive 
with a scale set by the radius of the S 3 . There is a single massless mode which is the zero 
mode of the temporal component of the gauge field. 

a = -L / A . (2.1) 

This mode is therefore strongly self-interacting even at arbitrarily weak 'tHooft coupling A. 
Consequently, one can, in the free theory, exactly integrate out all the massive modes and 
obtain an exact effective action for the mode a. This analysis was carried out in [8] and one 
obtains a unitary matrix model in terms of the holonomy 2 

U = e if3a , (2.2) 

where (3 = y is the radius of the thermal circle 3 . 

One finds that the gauge theory partition function (with U(N) gauge group and restrict- 
ing to adjoint matter fields) on S 3 x S l is given by, 



Z((3) = j [dU] 



cxp 



^^Tr(C/ n )Tr(f/^ 

.71=1 n 



(2.3) 



2 This unitary matrix model representation of the finite temperature partition function was given first by 
[6] based on enumeration of states in the free theory. See also [7] . 

3 The unitary matrix model for the free gauge theory was obtained earlier by Sundborg [6] by counting 
states of the free theory. See also [7] . 
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where the coefficients a n (T) are given, in terms of x = e 13 , by 

a n {T) = z B (x n ) + (-l) n+1 z F {x n ). (2.4) 

Here z B (x) and z F {x) are single particle partition functions of the bosonic and fermionic. modes 
respectively. They completely capture the field content of the gauge theory. The explicit 
expressions for z B (x) and z F (x), for fields of different spin are given in [8]. 

The above expression was derived at zero coupling where one has only a one loop con- 
tribution from integrating out all the massive modes. For weak 't Hooft coupling, one may 
continue to integrate out the massive modes and obtain a more general (and more complicated 
effective action for the holonomy U). The structure of the effective action is now [8] 

Z(J3,X) = J[dU]expS eff (U) (2.5) 

where 

1 k 

Seff(U) = J> K} (A,T,iV) M n TrC/ni ( 2 - 6 ) 

{«,} i=1 

with the integers rij obeying n « = an d the coefficients ct/ n .\, of a term with k traces, 
making their first appearance at (k — 1) loops in perturbation theory and consequently having 
a planar contribution starting with X h ~ 2 . 

Therefore, in perturbation theory, all the non-trivial low energy dynamics of the finite 
temperature theory on S 3 is captured by this unitary matrix model. It is the properties of 
this general class of unitary matrix models that we will study in this paper. 

However, one can make a further important simplification. The order parameter for the 
large N phase transition exhibited by these models (reviewed in the next subsection) is TrU. 
Consequently, one can also imagine integrating out all the TtU n (with n / ±1) and obtaining 
an effective action purely in terms of TrllTrW . This is not easy to carry out explicitly. 
Therefore one can consider toy models of the form [10] 

Z = j [dU]e N2s ^ (x) , x = -^TrC/TrC/t , ( 2 .7) 

where 

S eS (U) = ai(A, T)Tr[/TrC/t + bl ^P {TtUTyU^) 2 + (Trf/Trf/t) 3 + ■ • ■ (2.8) 

with S(x) being convex and S'(x) being concave. The simplest such model is the so-called 
(a, b) model [10] in which one keeps only the first two coefficients in the S eB given in Eq.2.8. 



Z(a u h) = J [dU] exp aiTrC/TrC/t + (TyUTxU^' 
where a\ and b\ are functions of temperature T and A. 



(2.9) 
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2.2 Eigenvalue Density Analysis at Large N 

The above unitary matrix models can be analysed using standard techniques in the large N 
limit. We briefly review the results [8] [10] in this subsection. 
One introduces the eigenvalue density 

1 N 

i=i 

where the holonomy matrix U takes the diagonal form 

U = diag(e i0i ). (2.11) 

Let us start with the free partition function Eq.2.3. It can be expressed in terms of a 
functional S[a(9)] 

Z{0) = J[Da]e N2s ^\ (2.12) 

where 

S[a(0)] = J d6 l j d6 2 a(6 1 )a(6 2 )V(e 1 -e 2 ). (2.13) 
Here, the two body potential V(9, T) is given by 

oo 

V(0,T) =ln[2] + ^-(l-a n (T))cos(n0). (2.14) 



n=l 



The general case for arbitrary a n is actually quite cumbersome to analyse. A self consis- 
tent method was given in [8]. See also [17] for a more general method. As mentioned above 
the crucial order parameter is Trf7, therefore we will often concentrate on the case where we 
keep only the terms with TrC7, TrU^ . In other words, we set a n = for n > 1. For this case 4 
where only oq ^ 0, one can explicitly obtain the saddle points for the above functional 2.13. 
One finds the following 

• For oq < 1, the minimum action is for the eigenvalue density 

a{6) = i-. (2.15) 

The free energy is zero for this configuration (to order A^ 2 ). 

• For aq = 1, there is a continuous family of minimum action configurations (labeled by 
a parameter £) for which the eigenvalue distribution is, 

a(6) = -!-(l + 2£cos<9) < 2£ < 1 . (2.16) 
All these configurations also have free energy zero. 



4 One can estimate that a n ~ for n > 1, at the phase transition temperature. Even at higher temperatures, 
the contributions of the higher a n is typically small. 
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• For a± > 1, there is a new saddle point whose eigenvalue distribution function is given 

by 



^^^MtH^i^i (2 - i7) 



where, 



^(yhW 1 -^^ (2 - 18) 



Note that this configuration is gapped unlike the above ones. The free energy for this 
configuration is (expressed in terms of £) 



-N 2 T 



e - i ln(2«) - i 



(2.19) 



Thus we see that there is a first order phase transition at a\ = 1, which corresponds to a 
temperature T = Th- 

This was for zero coupling. Perturbatively, we have more complicated matrix models 
Eq.2.6. Restricting to models of the form Eq.2.7, such as the (a, b) model, we can once again 
carry out a saddle point analysis of the eigenvalue density, using a Hartree-Fock approach. 
We now review the results (see [10] for more details) 5 . 

The large N saddlepoint equation for the model 2.7 is 

j<r(6)d<j)Cot (^~Y^j = 2Se//Kkisni 9 (2.20) 

where 

The solutions of this saddle equation are given by, 

*1 = Seff(^l>l , < (71 < \ (2.21) 



and 



2 



The phase structure is as follows: 

• (i) For sufficiently low temperature the only possible solution is of Eq.2.21 

o-i = . (2.23) 

This stable saddlepoint in fact has the uniform distribution Eq.2.15. It is identified 
with the thermal AdS saddlepoint on the gravity side. 



5 One can also study the general models by expressing them in terms of a suitable transform [9], [10], [17] 
of the one plaquette matrix model [15] [16]. This is particularly useful when studying the vicinity of points 
where the large iV expansion breaks down. 
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• (ii) At a higher temperature To, we find two new solutions, now of Eq.2.22 and thus 
with non-zero o\. The eigenvalue distribution for both these distributions are of the 
same form as Eq.2.17 (with different values of the parameter £). One of these is stable 
and the other unstable. On the gravity side, they can be identified with the small black 
hole (SBH) and the big black hole (BBH) respectively. 

• (iii) There exists a temperature T\ > To where the stable saddle points of (i) and (ii) 
exchange dominance. This temperature corresponds to the Hawking-Page temperature 
where the BBH has a lower free energy than thermal AdS in the semi-classical gravity 
path integral. 

• (iv) At some temperature T c > Tq, the eigenvalue distribution of the unstable saddle- 
point in (ii) changes from the gapped one in 2.17 to the ungapped one in 2.16. This 
Gross- Witten-Wadia(GWW) like phase transition has been identified by the authors of 
[10] with the black- hole string transition (see also [11]). 

• (v) And finally there exists a temperature Th, the Hagedorn temperature, when the 
unstable saddlepoint merges with the saddle point in (i). Above Th the saddlepoint in 
(i) becomes tachyonic. 

It is quite remarkable how the general class of matrix models 2.7 captures all the detailed 
qualitative features of the Hawking-Page phase diagram. This has been subsequently gener- 
alised to the case where one has a charge or chemical potential [21] [22] (see also [23]). One of 
the main new features [21] argued in the case of fixed charge is that we have a term in S e ff(U) 
of Eq.2.8 which is of the form ln((Tr?7Tr?7^). This has the effect that we no longer have the 
saddle point (i) above with uniform distribution. This matches with the known feature of the 
charged case that we never have thermal AdS as a saddlepoint. The only saddlepoints are of 
the form 2.17 and 2.16 and correspond to different kinds of black holes, small and big, stable 
and unstable. We refer the reader to [21] [22] for the details of the matching with the gravity 
phase diagram. 

3. Exact Solution at Finite N 

In this section we will obtain an exact expression for the gauge theory partition function 2.3 
of the free theory. As we will see the method of solution can be straightforwardly generalised 
to the general case described by Eq.2.6. 

Starting with the matrix model which captures the free gauge theory, 




(3.1) 



we can expand the exponential to obtain for the integrand 




(3.2) 
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Here, 

zS = U k ^ ki ( 3 - 3 ) 

j 

and 

oo 

T s (U) = I[(TrW) k '. (3.4) 

It is convenient to write this in terms of group characters using the Frobenius formula, 

T,:(U)=^2xR(C(k))Tr R U, (3.5) 

R 

where XR(C(k)) is the character of the conjugacy class C(k) of the permutation group 6 Sk, 
(K = ^2jkj), in the representation R of U(N). 

Now we can carry out the integral over the holonomy using the orthogonality relation 
between the characters of U(N) 7 , 

J [dU] Tr R (U) Tr R ,(tf) = $rr> ■ (3.6) 

Therefore we obtain 




This is an exact expression for any (5 and N for the partition function of the free gauge theory. 
In the next section we will further analyse the properties of this solution. However, it should 
be noted that the answer is completely explicit. The sum over representations of U(N) can 
be labelled by Young Tableaux with N rows and arbitrary numbers of boxes in each row. 
The characters of the conjugacy class are determined recursively by the Frobenius formula 
[18]. Explicit expressions for the most general case are not simple but have been given in the 
literature [19]. 

In the strict N = oo limit where the sum over representations is unrestricted, we can use 
the group theory identity for the orthogonality of characters of different conjugacy classes 
(see for instance, [18] pg.110) to obtain 

^09) = z; 11°? = IK 1 - ■ ( 3 - s ) 

k j j 

This agrees with the exact TV = oo answer derived in [8]. 

In the special case where only a\ / 0, the exact answer 3.7 is given by the simpler 
expression, 

oo 1 
fc=0 R 

6 Recall that a conjugacy class of the permutation group can be labeled by a partition k = (fci, fe, • • •)• ^ is 
an infinite dimensional vector with kj being the number of cycles of length j. 

7 The invariant Haar measure [dU] that appears here has been normalized such that J[dU] = 1. 
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In this case the only conjugacy class that contributes is the trivial or identity class consisting 
of k one cycles. The character of this class is nothing but the dimension dn(Sk) of the 
representation R for the permutation group Sk- 
it is clear that the above method of solution can be straightforwardly generalised to 
the matrix models 2.6 which describe the perturbative gauge theory at finite temperature. 
We can similarly expand the exponential and carry out the unitary integrals after using the 
Frobenius relations. The general answer for the finite N matrix model can once again be 
explicitly written though the actual expressions will now be more cumbersome. As a special 
case consider the (a, b) model 2.9. 



Z(a u h) = J [dU] exp cuTrC/TrC/ 1 " + (TvUTtU^ 



(3.10) 



Expanding the exponential as before and using a similar logic as before, we can write the 
partition function as, 

fc/2 k-2l h l r,, j2 (c \ 

k=0 1=0 y ' R 

For the more general case of models Eq.2.7 with 

oo 

fc=l 

We can again go through the same steps to obtain 

oo / k t \ 

*({<*}) = £ II ( jy2fc"/-i)fc I ) E ; K = J2lk, (3.13) 



R 1=1 



These are the cases whose large N limit we will be analysing in some detail in what follows. 

Finally, we should remark that matrix models of the form Eq.3.1 also appear in the 
counting of BPS states [31]. In fact, it is not difficult to use certain standard identities 
for the completeness of group characters to evaluate the answer Eq.3.7 in the case where 
a n = a i- One reproduces the usual generating function for the half BPS states given in [31]. 
It would be interesting to use the exact answer for finite N to evaluate some of the partition 
functions/indices of interest in N = 4 Super Yang-Mills theory [31]. 



4. Taking the large iV Limit 

We start by analysing the large N behaviour of the exact answer (for the free gauge theory) 
given in 3.7 8 . We should be able to see the large N phase transition that was obtained from 
the analysis of the eigenvalue density (reviewed in Sec. 2). We can see, in general, from the 

8 We will generalise to the interacting case in Sec. 4. 4 
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form of the solution that, in the large N limit, there is likely to be a dominant representation 
contributing in the sum over representations. Essentially, this can be viewed as a statistical 
mechanical system in which the group characters behave like an entropy contribution (roughly 
favouring representations with a large number of boxes in the Young Tableaux ). And the ctj 
are the Boltzmann suppression factors which disfavour representations with a large number 
of boxes. The balance between them leads to a dominant representation at any particular 
value of the temperature. A large N phase transition would occur when the nature of this 
dominant representation undergoes a qualitative change as one varies the temperature. As 
mentioned earlier, the Douglas-Kazakov [12] [14] phase transition and its generalisations in 2d 
Yang-Mills theory can be understood this way. We will now see all these features explicitly 
in our matrix models, specialising for simplicity to the special case where only a\ 7^ 0. As 
mentioned in Sec. 2.2, this case captures all the essential physics of the finite temperature 
theory. 

In the special case where a n = for n > 1, the exact answer is given by Eq.3.9 

00 ^ 

^(/3) = EE^[^)M • (4.1) 

k=0 R 

To proceed, we will write the sum over representations of U (N) in terms of the number of 
boxes of the corresponding Young Tableaux 

00 00 

E^E E <5( s *=i™i- fc ) with ni > n 2 > ... > n^v > , (4.2) 

R k=l {rii}=0 

where rij is the number of boxes in the i th row of Young Tableaux (there being only N rows 
for a representation of U(N)). Also k is the total number of boxes in the representation. 
Therefore the partition function reads as 

00 00 1 

= E E ft? ^ s ^ 2 «i s vlim - k ) ■ ( 4 - 3 ) 

fe=0{n l }=0 

The dimension d,R(Sk) is given by the Frobenius-Weyl formula [18] 

Kj 

where, 

hi = m + N -i , (4.5) 

with 

hi > h 2 > ... > h N > . (4.6) 
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4.1 The Continuum Limit and Saddlepoint Equations 

In the N — > oo limit we can define, following Douglas and Kazakov [12], continuous functions 
which describes each young tableaux 

n ( x ) = ^> K x ) = ^ , x = Jv 1 ^ 

where x £ [0, 1]. The function n{x) or equivalently h(x) captures the profile of the large N 
Young tableaux. In this limit Eq. 4.5 can be written as, 

h(x) = n(x) + 1 - x . (4.8) 

Note that the condition n\ > ri2 > ■■ > implies a strict monotonicity for h(x) 

h(x) > h(y) for y > x. (4.9) 

In this limit the total number of boxes in a Young Tableaux is given by, 

k = J2 n i N 2 (^J q dx [h(x) + x - l]j =N 2 ^jf dxh(x) - . (4.10) 

Since k' = dxh(x) — \ will generically be 0(1), we see that the number of boxes k = N 2 k' , 
in a generic representation, is of the order of iV 2 in this limit. 

The partition function Eq.4.3 can be written, using Eq. 4.4, as 

Z({3) = J^exp I ln(fc!)+lnai + 5^1n|^-/ij| -2^1n(^!) J . (4.11) 

{hi} \ i+j i ) 

In the large N limit, using Stirling's approximation for the factorials and Eq.4.10, the partition 
function can be expressed as 

Z = j [dh(x)]e- N2s ^, (4.12) 

where, 

- S eS = [ dx I dy In \h(x) - h(y)\ 
Jo Jo 

- 2 / dxh(x) In h(x) + k' In ( ai k') + k' + 1 . (4.13) 
Jo 

Recall that k! = j} dxh(x) — \. 

Now we are in a position to carry out a saddlepoint analysis for the effective action 
functional S eS [h(x)] (4.13). Varying S cB with respect to h(x), we obtain the saddlepoint 
equation, 

=lnh{x)--ln[ ai k'}. (4.14) 



f 

Jo 



h(x) - h(y) v ' 2 
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Introduce, again following [12], the density of boxes in the Young Tableaux u(h) defined 
By definition, it obeys the normalisation 



J hi 



dhu(h) = 1, (4.16) 



where the interval of support [hi,, hjj] of u(h) is specified by hi = h(l) and hjj = h(0). From 
the monotonicity of h(x) Eq.4.9, it follows that u(h) obeys the constraint 

u(h) < 1 . (4.17) 

In terms of the density u(h), the saddle-equation (4.14) can be written in the more 
familiar form, 



(4.18) 



where £ 2 = a\k' . Note that the parameter £ involves k' given by 

k' = dxh{x) ~\ = J hu{h)dh - ^ (4.19) 

which in turns depends on the (first moment of the) density u(h). We will therefore have to 
solve the equation self-consistently. 

4.2 The Saddlepoint Densities 

The solution to the integral equation 4.18 for the Young Tableaux density u(h) is obtained 
along similar lines to the usual solution for the eigenvalue density. The main point to addi- 
tionally take into account is the presence of the constraint Eq.4.17. Thus we will find that 
the solutions to 4.18 are of two different kinds depending on the value of the parameter £. 

Solution Class 1: 

0<u(h)<l; h£[q,p}. (4.20) 

A typical representation corresponding to such a Young Tableaux density is plotted in figl. 
There are always a nonzero number of boxes in each row. 

Solution Class 2: 

u(h) =1 h £ [0, q] 

= u{h) h G [q,p] (4.21) 

with < u(h) < 1. A typical young tableau for a solution of this class has been plotted in 
Fig.l. The representations are such that a finite fraction of the rows are empty. 
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n l (x=0) 



n N (x=1) 




Solution Class 1 



Solution Class 2 
Figure 1: Young Tableaux 



As we vary £, the constraint 4.17 will come into play and one will have to switch from 
one of the branches to the other. At this point, as we will see explicitly, there will be a 
non-analyticity, for example in the free energy and we will have a large N phase transition. 

We will now solve the saddle-equation in the conventional way 9 by introducing the resol- 
vent H(h) defined by 

„ <h>) 



H(h) = \ 

Jh L 



dti 



h-b! 



(4.22) 



The resolvent has the following properties, 



• (i) It is analytic in the complex h plane with a branch cut along the positive real 
interval (hi-, hjj). 

• (ii) It is real for real positive h outside the interval. 

• (iii) H(h) ~ ^ + (k' + (as h — > oo). This follows from the moment expansion of 
the resolvent at large h and using Eq.4.19. 



(iv) H(h + ie)+H(h-ie) =2 In 



for real h . 



• (v) u(h) = - J;, [H(h + ie) - H{h - ie)] for h G [q,p]. 

One can therefore solve for H(h) in terms of its real part by writing it as a contour 
integral. In fact, the equations we need to solve are very close to the equations that arise in a 
class of matrix models studied by Kazakov, Staudacher and Wynter[13] (see also [25] for the 
solution of a similar equation). We now exhibit this solution in both the classes mentioned 
above. 



9 For a recent review see [20]. 
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Figure 2: Contour for Solution class 1 



4.2.1 Solution Class 1: 

In branch 1 using the ansatz 4.20 the saddle equation becomes, 

~K 



In 



he[q ,p]. 



(4.23) 



The resolvent (Eq. 4.22) in this branch is given by (see [13], [20] for a general discussion), 

ds In (a/f) 



H(h) = -y/(h-p) {h-q)j 



(4.24) 



2vri {s - h )^ {s - p ) {s - q ) ■ 

The contour of integration is shown in Fig.2. Carrying out the contour integration we obtain 
the resolvent for class 1, 

"2 h 2 - (y/p - ^qf h + 2 q p - 2 (h + ^qp) ^{h - p)(h - q) 



H(h) = In 



(4.25) 



Using this, we can readily find the discontinuity and thus the Young Tableaux density 
function u(h), 

>2" 



u(h) = —cos 1 

7T 



—cos 

7T 



2^ ' 2£h 

h + Z-l/2 



for h £ [q,p] 



2^h 



(4.26) 



Fig.3 shows the plot of u(h) vs. h for this solution class. 

The support of u(h) as well as k' is determined by expanding H(h) for large h and 
matching with property (iii) of H(h) listed above. 

'(v^ + v^) 2 l , (v^-v^) 21 



H(h — ► oo) ~ In 



4£ 



+ 



T + (^4U. (4.27) 
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Figure 3: Plot of u(h) vs. h for Solution Class 1. The value of £ increases from 0.5 as one goes from 
the leftmost graph to the right. 



We therefore obtain 



and 

which implies 



k' = y/qp + - 



(4.28) 
(4.29) 

(4.30) 

4(-l <431) 
using the definition £ 2 = a\k' . Since q is a real positive quantity, Eq.4.28 implies that this 
solution branch exists for £ > \. From Eq.4.31, we therefore conclude that this class of 
solutions only exist for oq > 1. 

4.2.2 Solution Class 2: 

Using the ansatz 4.21 the saddle equation becomes, 



ai 



4e 



In 



h-h' 

The full resolvent H{h) now takes the form, 

//!//! hi 



In 



h 



h 



h 



+ 



dh 



where, h G [q,p] . 
h-h' ' 



H(h) can once again be written as a contour integral, 



H{h) = In 



h 



h — q 



-V(h-P) (h-q)j 



ds ln(s/£) -\n[s/{s-q)} 



2lTi ( S - h )y/(s-p) (i 



(4.32) 



(4.33) 



(4.34) 
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q p 



Figure 4: Contour for Solution Class 2 





Figure 5: Plot of u(h) vs. h for Solution Class 2. The value of £ decrease from 0.5 to as one goes 
from the rightmost to the leftmost graph. 



The contour is shown in fig 4. Carrying out the integration gives the answer 
H(h) = In 



In 


'h~ 


+ In 




+ In 




A. 







h + ^qp- V(h-l) 2 -4e 



h 



qp + y/(h - l) 2 - 4^ 2 



Hence, the Young Tableaux density is given by, 

u(h) = icos -1 

7T 



h-1 



(4.35) 



(4.36) 



Fig. 5 shows the plot of u(h) vs. h for this solution class. 

As before, expanding H(h) for large h we find the values of q, p and k' as follows, 

q = 1 - 2£ , (4.37) 
p = 1 + 2£ (4.38) 
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and 

k' = • (4.39) 

From the definition a\k! = £ 2 we obtain 

either £ = 

or,cn = 1. (4.40) 

The former implies the uniform distribution 

u(h) = l he [0,1]. (4.41) 

This is therefore a saddlepoint for any value of a±. This is in fact the density corresponding 
to the trivial representation rij = 0. 

The latter corresponds to a family of saddlepoints labelled by £ which exists only at 
ai = 1. From Eq. (4.37) it is clear that this family exists for £ < 1/2. 

To summarise, we see that there are three different saddlepoint configurations of Young 
Tableau densities. 

• The trivial representation corresponding to the uniform distribution for u(h), Eq. 4.41. 
This exists for any value of a±, i.e. for any temperature. 

• The continuous family Eq. 4.36 which exists only for a\ = 1. Here a finite fraction of 
the rows of the Young Tableau are empty. 

• The representation Eq. 4.26 which exists only for a\ > 1, i.e. only for high enough 
temperature. Now all the N rows of the Young Tableaux are filled. 

We see that the saddlepoints are exactly in correspondence with the saddlepoints of the 
eigenvalue density, as summarised in Sec. 2. 2. We will now indeed verify that the large N free 
energy of these saddle points is also exactly the same as that seen from the eigenvalue density 
analysis. 

4.3 The Free Energy of the Free Theory 

In the large N limit the free energy of the partition function Eq. 4.12 is given by, 

F = —T\nZ = N 2 TS° a . (4.42) 

where S® s is the value of effective action at the (dominant) saddlepoint. The effective action 
given in 4.13 can be re-expressed as a functional of u{h) 

-S eS = dhl dh'u{h)u{h')\n\h-h'\-2 dh u{h) h \n[h]+k' + l+k' ln[ ai k'}. (4.43) 

We only need to evaluate this functional on the different saddlepoint configurations we have 
found. In evaluating the expressions it is useful to use the corresponding saddlepoint equations 
to eliminate the quadratic term in u(h) in the effective action. We thus obtain: 
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Solution Class 1: Using the saddle equation Eq. (4.23), S e3 becomes 

S% = - £dhu(h)h\nh+ ^-^lnK] + i+d , (4.44) 
where the constant C\ is given by, 

Cl = l -Hi/2\+i- 1 -. (4.45) 

Evaluating the integral (details can be found in Appendix A) gives finally for the effective 
action, 

1 , , 1" 



F = N 2 TS°„ = —N 2 T 



< 0. (4.46) 



This exactly agrees with the free energy computed in [8] for the deconfined phase sad- 
dlepoint which was quoted in Sec. 2. 2. 

• Solution class 2 : 

Here the free energy is easy to compute since we have a continuous family labelled by 
£ which are all saddlepoints and thus must have the same free energy when a± = 1. 
In particular the constant configuration Eq.4.41 is a limiting member of this family for 
which the effective action is readily computed to be zero. Therefore the entire family 
of saddlepoints in Eq.4.36 must have zero free energy. This is explicitly verified in 
appendix A. Once again this matches with the results of the eigenvalue analysis for the 
confined phase saddlepoint. 

We see from this that there is an exchange of dominance of the saddle points at a\ > 1, 
where one gets a new saddlepoint Eq.4.26, which has less free energy compared to the uniform 
density saddlepoint that exists for all values of a\. Therefore as the temperature increases we 
get a transition when the saddlepoint switches giving rise to the confinement-deconfinement 
transition. We see that this is nothing bu the analogue of the Douglas-Kazakov transition in 
our approach. We have thus reproduced the usual saddlepoints as well as the phase diagram 
of the zero coupling Yang-Mills theory from our method of taking the large N limit of the 
exact solution. 

4.4 Extension to Non-zero Coupling 

The general interacting unitary matrix model Eq.2.6 for perturbative gauge theory can also 
be solved exactly using the technique of Sec. 3 and a large N limit can then be taken along 
the lines described in this section. However, the analysis is going to be technically much more 
involved. But as described in Sec. 2, the essentials of the physics is in any case captured by 
models involving only TrU. In particular, the (a, b) model Eq.2.9 already does a good job in 
getting the detailed form of the Hawking-Page phase diagram[10]. Here we will sketch how 
its exact solution Eq.3.11 shows all the features described in Sec. 2. 2 when we take the large 
N limit as described in this section. 



- 19 - 



Let's write the answer Eq.3.11 as 



Z ( ai , h)=j2i: N 2 ln{k _ 2l)l e ~n^ 5{k - ^ 

k=o 1=0 v ' {m} 

00 d 2 (Su) 

= E ^ E -^j^^ - ^ni) , (4.47) 
k=o {m} 

where, as usual, {raj} label the number of boxes in the Young Tableaux and 

Since the total number of boxes k is of order iV 2 , we can replace the sum in Eq. 4.48 by the 
saddlepoint value. Doing this gives 

k 



k ( a\ \ _kx _ k„ k 



where x is determined by the equation 

x 26i k' 

(1 - xf ~ ~af 

with k = N 2 k' . Therefore the partition function 



(4.50) 



{rii} 

(with k = rii) takes essentially the same form as Eq.4.3 except that we have a\ instead of 
a-y. 

We can now take the large N limit as before to obtain 

Z{a x M) = J [dh(x)] e-^ff, (4.52) 

where S e3 is the same as in Eq.4.13 with the replacement of a\ by a\. 

Since d\ depends on k' or h(x) (see 4.49 and below), the saddle equations are modified a 
bit. Taking into account this additional dependence on h(x), the saddle equation becomes, 



Jh, 



h(x) 

h-h! ~ [~J~ 



where £ 2 = ai£/e2. (4.53) 



Since the saddle equations are of the same form as Eq.4.18, (with the replacement of £ by 
£) the saddlepoint configurations for the Young Tableaux density are also the same in form. 



3 The extra multiplicative factor of -| plays only a subleading role in the large N limit. 
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Namely, we obtain the three different configurations of Sec. 4. 2. In fact, we can redo this 
analysis for the class of models Eq.2.7 for which the exact solution was given in Eq.3.13. This 
is performed in Appendix B. We see from the analysis there that the saddlepoint equations 
give once again the same saddlepoint configurations. Moreover, we obtain the same Hartree- 
Fock equations as Eq. 2.21 and Eq. 2.22. Thus the phase diagram turns out to be the same 
as that given by the eigenvalue density analysis. For instance in the case of the (a, b) model 
we find 

• A low temperature saddlepoint which is characterised by £ = which is the uniform 
distribution corresponding to thermal AdS. This has zero free energy. 

• Then there is a saddlepoint of the form Eq. 4.36 when £ 2 = k' i.e. a\e^ = 1 and £ < ^. 
This implies that 

~e = ^ < - . (4.54) 

^ 26i - 4 y J 

This is actually the unstable saddlepoint corresponding to the small black hole (in the 
phase where it is to be viewed as an excited string state) and has positive free energy 

F = AT 2 r (1 ~ ai)2 . (4.55) 
46i 

This saddlepoint exists in a temperature range T c < T < Th- 

• Finally there is the saddlepoint of the form Eq. 4.26 which obeys 

ai _ 4g 
1 - x 4| - 1 

and k' = f - - . (4.56) 

The two solutions to this equation give rise to the BBH as well as the SBH (in the 
actual black hole regime). The BBH solution exists for all temperatures greater than 
the minimum To for which this solution exists. While the (gapped) SBH solution exists 
in the interval T Q < T < T c . At T c which corresponds to £ = ^, this solution goes over 
into the ungapped solution 11 . The free energy in this phase is given by, 



-N 2 T 



1 . , r, 1 , ( 1 N 4 



(4.57) 



One of the points to note in our analysis is that the saddle configurations for u{h) are 
all of the same form as in the free theory. It is only that £ is a different function of the 



11 Therefore the GWW transition identified in [10] with the black hole-string transition is the same as the 
Douglas-Kazakov (DK) transition in our approach. A similar thing was seen in [14] where the DK transition 
of 2d Yang-Mills was mapped onto a GWW-like gapped to ungapped transition in terms of the eigenvalues of 
Wilson loops. 
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temperature. This is a mirror of the same phenomenon in the eigenvalue density analysis 
that the functional form of the saddle configurations of a(9) are not changed as one turns 
on the perturbative coupling. This robust character of the saddle point configurations is a 
positive indication in trying to extract universal features from these results. 

Finally, we should mention that the extension to non-zero charges also follows in a 
straightforward way. As argued in [21], the effective matrix model has a logarithmic term 
which results in there no longer being a uniform eigenvalue density saddle point corresponding 
to thermal AdS. In our approach, as argued in Appendix B for a general matrix model of the 
form Eq.2.6, we find the same saddle point equations as in the eigenvalue analysis and thus 
the same phase diagram. 



5. Free Fermionic Phase Space Description 



5.1 Relation between the Young Tableaux and Eigenvalue distributions 

Our analysis of the exact answer has been very different from the usual eigenvalue analysis 
reviewed in Sec. 2. It turns out, rather remarkably, that there is nevertheless, a simple rela- 
tionship between the saddlepoint configurations u(h), in both the high and low temperature 
phases, with the corresponding saddlepoint eigenvalue densities. We will now describe this 
relation. 

Consider first the low temperature saddlepoint u(h) = 1. The corresponding saddlepoint 
for the eigenvalue density is a(6) = From the graph of these two distributions, we notice 
that they are functional inverses of each other (flip the horizontal and vertical axes of one to 
get the other). In other words, we can make the identification 



6 

u = — 

TT 



(5.1) 



This case may seem a little trivial, so let us consider the family of saddlepoints that 
correspond to the unstable small black hole Eq.4.36 



u(h) = — cos" 

TT 



1 



2£ 



h £ [q,p]; 2£ < 1 . 



(5.2) 



together with u(h) = 1 for h £ [0,p]. Applying relations 5.1 to this case, we get immediately 
a(0) = ±(1 + 2£cos<9) which is the same as Eq.2.16. Thus the identification holds in this 
case. The two distributions are again functional inverses of each other. 12 

We finally come to the non-trivial saddlepoint corresponding to the big black hole Eq.4.26 

h€[q,p]. (5.3) 



u(h) = —cos 

TT 



h 



2£ 



12 A similar relation was also found in the case of 2d Yang-Mills theory on the cylinder[14] though a phase 
space interpretation was not made. 
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In this case, we have to be more careful. We see from the plot of u(h) that the functional 
inverse is ambiguous. For a given value of u, there are two values of h. This can be directly 
seen from the fact that Eq.5.3 implies the quadratic relation 

h 2 -[l + 2ZcoB(iru(h))]h+(z-^ =0. (5.4) 
If we take the difference between the two solutions h + and h- , we obtain 



h+-h- = 2^2^/l -2£sin 2 (^) cos (™) . (5.5) 
We see that if we define 



,2#o 1 



2 = 2? (5 ' 6) 



and modify the identifications Eq.5.1 to 



u = — 

IT 

h " = a{0) , (5.7) 



2vr 

then we obtain precisely the eigenvalue distribution in Eq.2.17 



a(9) = Jsin 2 ^ - sin 2 - cos [ -) . (5.8) 




5.2 Fermionic Phase Space 

The relations Eq.5.1 and 5.7 between the saddlepoint eigenvalue densities and the young 
tableau densities have a very natural interpretation in terms of a free fermionic picture. This 
fermionic picture is suggested by the fact that the eigenvalues of the holonomy matrix behave 
like fermions. At the same time, the representations of U(N) also have an interpretation in 
the language of non-interacting fermions with the number of boxes of the Young tableaux 
being like the momentum (See [24], for example). This suggests that the eigenvalue density 
is like a position distribution while the Young Tableau density is like a momentum distribu- 
tion 13 . Therefore, it is natural to consider a phase space distribution which gives rise to these 
individual distributions. In the classical (i.e. large N) limit, we can describe this system of N 
fermions in terms of an incompressible fluid occupying a region of the two dimensional phase 
space (see [30] for a recent overview and references to the large literature on the subject). 

Therefore let us assume that the saddlepoints are all described by some configuration 
in phase space, i.e. some region R of the two dimensional plane such that the phase space 
density p(h,9) obeys 

p(h,9) = ±-; (h,0)€R 

= ; otherwise. (5.9) 

13 In fact, u{h) can be viewed as a plot of the fermi distribution of momenta. The uniform density saddlepoint, 
for instance, corresponds to a fully filled fermi sphere. 
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Figure 6: Phase space distribution for Thermal AdS: £ =0 

We can then define the partial densities 14 

u(h) = ! p{h,9)d6 

J —71 
/"CO 

<t{0) = / p(h,0)dh (5.10) 
Jo 

where the first integral is at constant h and the second at constant 9. Note that 

J P (h,9)dhd9 = 1 . (5.11) 

Let us now take the boundary of the region R to be defined by a curve 15 C(h, 9) = O.We 
now see that there can be different situations depending on whether the solution h = h{6) is 
single valued or multiple valued. Thus for a single valued h = h(9), we see from Eqs.5.9 and 
5.10 that the identification Eq.5.1 follows. When we have multiple values for the solution h{9) 
then the identification between the different densities is a little more non-trivial. For instance, 
when we have two solutions h+(d) and h-(0) with h+ > h-, then it follows from Eqs.5.9, 5.10 
that the relation between the young tableaux density and the eigenvalue densities is that in 
Eq.5.7. Thus we have an interpretation for the relationships that we observed in the previous 
subsection. 

14 The measure factor that appears in Eq. 5.10 and Eq. 5.11 suggests that h is related to the usual polar 

2 

coordinate r by h = ^- . This redefinition is quite natural from the point of view of free fermionic phase space 
where a similar change of variables is made. See [30] and [3]. 

15 The symmetry of the effective action under U — > implies the region is symmetric under 9 — > —9. 
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Figure 7: Phase space distribution for SBH in ungapped phase: £ =0.4 



Figure 8: Phase space distribution at GWW transition point : £ = 0.5 

We therefore see that the large ./V saddlepoints of the gauge theory effective action, 
which correspond to the Thermal AdS, the small black hole and the big black hole can each 
be thought of in terms of a particular configuration in a free fermionic phase space. There is 
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Figure 9: Phase space distribution for SBH (in the gapped phase)/BBH 



a particular shape associated to each of them. This shape, which is determined by the curve 
C(h,0) = 0, can, in general, only be inferred from the knowledge of both o~(9) and u(h). For 
the different saddlepoints that we have discussed, the particular shapes are given in Figs. 6- 
9. We have plotted the regions in polar coordinates (r, 9) after making the identification in 
footnote 13. 

The shape corresponding to thermal AdS is a disk h{6) = 1 (Fig. 6) which is exactly 
like that in the Lin, Lunin, Maldacena description for global AdS. For the big black hole we 
have a shape Fig. 9 which corresponds to a double valued h(9). The equation of the curve 
is given by, Eq. 5.4 with iru(h) replaced by 9. Note that the origin of phase space is not 
contained in this region. For the unstable saddlepoint (SBH) we have the black hole regime 
where the shape is qualitatively the same as the BBH, except that it is closer to the origin. 
At the temperature T c , the shape continuously changes to that in Fig.8. The excited string 
state beyond this transition occupies a region which includes the origin given by the curve 
h{9) = l + 2£ cos 9 (Fig. 7). We observe that the shapes are qualitatively of two different kinds 

2 

with the limiting shape h = ^- = (1 + cos 9) at the GWW transition point that separates the 
two classes. Note that the boundary of the limiting configuration is a separatrix between two 
different kinds of trajectories in the fermi sea 16 . 



'We thank S. Wadia for useful discussions on this point. 
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6. Conclusions 

We have studied thermal gauge theory by evaluating its partition function at finite N. Taking 
the large N limit of the full answer gives us a new perspective on some already known facts 
about the phase diagram of the gauge theory. The various phases and the transitions between 
them can be viewed in terms of different dominant representations of U(N) that contribute to 
the partition function. They are obtained as saddlepoints of an effective functional, S[u(h)] 
(like that given in Eq.4.43) of a young tableau density. These representations can be viewed as 
excitations of a one dimensional fermi sea whose ground state corresponds to (thermal) AdS. 
In fact, a complete picture emerges when we combine the saddlepoint representation u(h) with 
the saddlepoint a{9) of the eigenvalue density analysis. One obtains a phase space picture 
which neatly combines both these answers into a single two dimensional region with uniform 
density p(h,9). An important point to emphasise is that in general, one needs both u(h) and 
a (9) to obtain the curve C(h,9) = which delineates the region R of phase space which is 
occupied. In a way, as we have seen, u(h) is more basic since it contains the information to 
reconstruct cr(9) through Eq.5.7 but not vice versa. 

This suggests that it is more natural to look for a description of the dynamics directly in 
terms of the phase space density p(h, 9). In particular, it would be nice to find a phase space 
functional S[p(h, 9)} which reduces to the two different functional S[u(h)] and to S[a(9)} on 
appropriately integrating out. It would seem that the formalism reviewed for example in [30], 
would be the appropriate one. 

Such a formulation is likely to help in moving towards the goal of reconstructing the local 
theory in the bulk with all its redundancies. Note that the configurations we are considering 
in the thermal history are all S" 3 invariant. Thus the only non-trivial directions in the bulk 
are those of the thermal S 1 and the radial direction. It is tempting to identify these two 
directions with the phase space directions of the fermions. In fact, the eigenvalues of the 
fermions can be viewed as positions on the T-dual to the thermal circle and the momenta 
should correspond to the radial direction as per the usual AdS/CFT correspondence. As a 
first step it would be nice to see how the topology of the bulk is encoded in the geometry 
of the phase space regions corresponding to the different saddlepoints. Note that we argued 
that the shape of the regions is fairly robust against coupling effects and so should be telling 
us something generic about the bulk geometries even when they are subject to all kinds of a 
corrections. 

It may seem a little confusing to talk about phase space without talking about time. But 
since we are describing euclidean gravity configurations in terms of euclidean gauge theory, 
time does not enter directly. A closely similar situation arises in the study of Euclidean 2d 
Yang-Mills theory on the torus. Its free fermion representation was interpreted [26] in terms of 
contributions from different "baby universe" saddle points of the Euclidean gravity partition 
function of the dual geometry (in this case a 4d extremal black hole in an N = 2 theory) [27]. 

As mentioned in the introduction, the LLM scenario was one of the inspirations for many 
of these ideas. It would be nice to see if there was a concrete connection to be made between 
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those cases and those studied here. At first sight they seem to be very different. However 
one possible link between them is that unitary matrix models very similar to those studied 
here give the partition functions that count BPS states. In fact, it would be nice to apply 
our finite ./V answers for these models to these countings 17 . 

Another interesting direction to generalise might be to the gauge theory on spatial sections 
which are for instance, S 1 x R 2 rather than S 3 . We now have a matrix quantum mechanics 
[28]. Toy matrix models for this case have been studied for example in [29]. In these cases a 
free fermionic description has been employed to study the matrix quantum mechanics. 
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Appendix 



A. Details of the Evaluation of the Free Energy 



A.l Free Energy for solution class 1 

The effective action is given by 

"h TJ rhu 



I dhj- 

hr J hf 



dh'u{h)u(h')ln\h-h'\-2 dh u(h) h\n(h)+k' + l + k'ln( ai k') (A.l) 

Jh L 



For solution class 1, hi, = q and hjj = p. u{h) has support only for q < h < p. 

Using the following saddle equation of motion we can replace the double integration by 
a single integration, 

l P dti \n[h - ti] u{h') = h\n l! 

J Q 



h + d 



Hence the effective action becomes, 

S° cS = ~ dhhu{h) \u[h] + i ln[£] - ln[£] + i + d 
The constant C\ can be evaluated from the equation A. 2 at h = p. 



+ 



^ dh'\n\p- h'} u(ti) . 

J a 



C\ =p + pln 
After some algebra we find, 

J" dhhu(h) In [h] = £ ln[£] + \ + \ 



and 



Finally we get, 



Cl = Iln[£/2]+£-i 



5 e ° ff = ln[£] - Ci 



A. 2 Free Energy for solution class 2 

The effective action for the solutions in class 2 is given by, 



-S cff = J P dh-f P dh' u{h)u{h')ln\h-h'\-2 [* ' dhu(h)hln(h) 
Jo J o Jo 

+ k' + 1 + k'\n( ai k'). 



(A.2) 



(A.3) 



(A.4) 

(A.5) 
(A.6) 



(A.7) 



(A.8) 
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Breaking the h integration into two pieces, we can write the effective action as, 

n ' dh'u{h)u{h') In \h - h'\ + [ P dh-f 1 dtiu{h)u(ti) In \h - h'\ 

Jq J 

fP 

2 dhu(h)h\n(h) + k' + 1 + k' \n( ai k'). (A.9) 



ri 

-S cB = I dh 
o 







In the second term on the right hand side, one can use the saddle equation Eq. (4.18) with 
hi = and h\j = p. This gives 

[ P dh-f P dh'u{h)u{h') In \h - h'\ = [* dhu(h) (h\n[h] - hhife] - h - C 2 ) , (A.10) 

Jq J Jq 

where the constant C2 is given by, 

C2 = p — p In 
After some algebra we get, 



p 



fdti u(h') \n\h-ti\ . (A.ll) 
J 

C 2 = ln[C] . (A.12) 



Therefore 

rp 



dh 

1 



;-f P dh'u{h)u{h')ln\h- h'\ = f dhhu(h)ln[h] 

Jo Jq 

- (1 + Mfl) (V + \ - y) + Ca(l - g). (A.13) 
Calculating the first term on the right hand side of Eq. (A.9) we get, 

f 1 dhl P dh'u{h)u(ti)\n\h-ti\ = -^ + q 2 ln[q] + f P dh!h'u{h') \a[h'] 
Jo J 2 J q 

- q+ T dh'u{h'){q - h') \n[ti - q] . (A.14) 

Jq 

Calculating other terms in on the right hand side of Eq. (A.9) we finally see that the on-shell 
effective action in second branch vanishes. Hence the free energy in this branch is zero. 

B. A Class of General Matrix Model Actions 

In this appendix we will generalize the result of section 4 to a generic effective action which 
is a function of Trt/Trf/t . We will expand the effective action in a power series in x, where 
x = TrC/TrC/t. 



jy2(n-l) jV"2n " 

n=l n=l 



Using this form of S aH , we can write 



m 00 



i ki=0 



,-Seff = TT V 1 i} (B 2) 
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The partition function is given by, 



Z = j [dU]e- s °«, 

can be evaluated using the methods of Sec. 3 and Sec. 4 

R 

m 

N ;=1 



(B.3) 



OO ki k 

Gti * 'Q>rr> 



{fei}=0 



fei! • -k m \ 



(B.4) 



z=i 



Define, 



oo ki k„ 
fti * 'CLt-, 



fc! 



{fej=0 1 m AT i=i 

oo 

e <3({fei}) j 

{fci}=0 



(B.5) 



where G({ki}) is given by, 

m m 

G{{ki}) =J2 k i m t a *] " Yj ( ki Mk] ~k t ) + k \n[k] W* ~ 1) M^]- ( B - 6 ) 



i=i 



i=i 



In the large N limit the partition function receives its dominant contribution only from the 
extremum value of G. Minimising G with respect to ki, remembering to introduce a Lagrange 
multiplier for the constraint ^ ifcj = /c, we get, 



i=i 



k; — flj/3 , 



(B.7) 



where k\ = jfe and j3 = Q~( a + 2lnN ) where a is the Lagrange multiplier enforcing the con- 
straint. 

(3 is determined through the following relation, 

oo 

1=1 



ukF 



i=i 



where <S"(/3) is given by, 



1=1 



Once we fix the undetermined multiplier, then G can be written as, 

-^2=k' ln[k'} -k'-k' ln[(3] + S(P) . 



(B.8) 
(B.9) 

(B.10) 
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Hence f{k) is given (upto multiplicative factors which are unimportant in the large N limit). 

/Wae^^-V^Cfl • (B.ll) 
So the partition function can be written as, 



a* 



^EfclEto, (B.12) 

k=0 ' R 

where a k = e khl i k> VP]- k + N2s (P) . Like before we will write the partition function in the 
following form, 

Z = J[dh(x)]exp[-N 2 S ca ], (B.13) 

where S eS is given by, 

-S eS = f dx dy\n\h(x) - h(y)\ 
Jo Jo 

- 2 f dxh{x)\nh(x) + k' In (a k') + k' + 1 . (B.14) 
Hence the saddlepoint equation is given by, 

where ^ is given by, 

= k'S'(p) . (B.16) 
For this generic effective action the solutions of the saddle equation are given by, 

(3S'((3) = (3 for P€[0,±] 

^-w^ for/?G[ i' 1] - (R17) 

For G [0, |], £ = VP and for G [|, 1], £ = 

These equations are identical to those obtained by the Hartree-Fock analysis of eigenvalue 
density Eq. 2.21 and Eq. 2.22. 
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